Local field factors in a polarized two-dimensional electron gas 



(N 
O 

o 

> 
O 

o 



i 

CO 

-4-5 
c3 



I 

C 

O 

o 



> 

(N 



(N 
O 

c3 



o 
o 



X 



Juana Moreno and D. C. Marinescu 
Department of Physics, Clemson University, Clemson, SC 
(Dated: February 1, 2008) 

We derive approximate expressions for the static local field factors of a spin polarized two- 
dimensional electron gas which smoothly interpolate between their small- and large-wavevector 
asymptotic limits. For the unpolarized electron gas, the proposed analytical expressions reproduce 
recent diffusion Monte Carlo data. We find that the degree of spin polarization produces important 
modifications to the local factors of the minority spins, while the local field functions of the majority 
spins are less affected. 
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I. INTRODUCTION 

In most of the many body theories of electronic sys- 
tems, the spin polarization in a magnetic field has been 
considered as a small parameter. This picture is cer- 
tainly true when the Zeeman splitting is much smaller 
than other relevant energies in the problem. However, re- 
cently, attention has been focused on materials where the 
Zeeman splitting dominates the energy spectrum, such as 
the diluted magnetic semiconductors. The large Zeeman 
splitting is due to the strong exchange interactions be- 
tween the itinerant carriers and the magnetic ions that 
generate a value of the effective gyromagnetic factor, 7*, 
up to hundreds of times its band value. As a conse- 
quence, even in weak magnetic fields these systems can 
be fully polarized.El Due to their easy polarizability these 
materials are promising candidates to accomplish spin 
dependent conduction in solid state devices ;□ hence the 
interest in developing an appropriate microscopic model 
for them. 

A realistic picture of the spin polarized electron gas 
hinges on finding an appropriate description of the many- 
body interaction, which has to incorporate the explicit 
spin dependence of the short range Coulomb repulsion. 
In this paper, we examine these short range effects in 
a two dimensional polarized electron gas. The self- 
consistent treatment of the exchange and correlation ef- 
fects has proven to be very important in understanding 
the physics of normal metals ,u but to our knowledge it 
has not been fully analyzed in spin-polarized systems. In 
addition, the relevance of the exchange and correlation 
effects increases as the dimensionality of the electron gas 
is lowered. 

We model the exchange (x) and correlation (c) hole 
around each electron by using sain dependent local field 
correction functions, G^' c (q, lu)u They describe the dif- 
ference between the true particle density and its mean 
field (RPA) counterpart. The spin dependence of these 
factors is a consequence of the fact that the microscopic 
interaction is different for up and down spin electrons. 
The microscopic origin of the local field corrections was 
elucidated by Kukkonen and OverhauserJj In a direct 
generalization of this approach, the effective potential V a 



experienced by an electron with spin a in the presence 
of external electric fl ||^(q, w)] and magnetic fields [B(q, u>)] 
can be written as:l 



V ff (q,w) 



(q,w) +7V • B(q,u)- 



>(q){[i 



Gt (q, w)]An(q, w) - a ■ Am(q, u)G a (q, w)\ 

(1) 

where An = Any + Any and Am are, respectively, the 
density and the magnetization fluctuation induced by the 
external fields, a are the Pauli spin matrices and v(q) is 
the bare Coulomb interaction, 27re 2 /q. The local field 
factor G+(q, u>) is the sum of parallel- and antiparallcl- 
spin effects, while G~(q, u>) is the difference of parallel- 
and antiparallel-spin effects: G^(q, ui) = G x a + G c aa ± 
G c - I 

W (7(7 

In a linear approximation, the density fluctuations are 
proportional to the effective potentials: Any = HyyVy 
(An j, = II|jVx)Ji'i where the proportionality coefficients 
are the polarization functions of the fully interacting elec- 
tron system. Since the interacting polarization function 
is generally unknown, an additional local field factor (G") 
is needed to relate the interacting polarization function, 
n CT(T , with the non-interacting one, 11°^., as inQ 
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I + 2 V (q)G«(q,^)II° ff (q, W ) 



(2) 



This parameterization of the modified polarization 
function leads to a renormalized expression of the lo- 
cal fields that determine the response functions, such 
that the complete local field functions are given as: 

/~r± /"<± 1 rm Q 

Quantitative calculations of many physical properties 
require the precise knowledge of the local field correc- 
tion functions. The determination of the frequency and 
wave vector dependence of the local field corrections is 
a very difficult problem which remains unsolved even in 
the case of the unpolarized electron system. Fortunately, 
the asymptotic values of the local field, factors can be 
obtained exactly in some limiting cases. ElE3Eil Numerical 
estimates of the response functions of the two and three 
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dimensional unpolarized electron gas have shown that lo- 
cal field factors smoothly interpolate between the asymp- 
totic small and large wave- vector behavior .L3 This feature 
is expected to exist also in the case of a spin polarized 
system, and, consequently, we use the asymptotic limits 
of the static local factors for large and small wavevector 
as a starting point in deriving their approximate expres- 
sions across the whole spectrum of momentum. 

The fundamental parameters of the problem are the 

coupling strength, r s = a* B / ^ ir{n^ + n{) = a%/\f~Kn$3 

and the spin polarization, £ = (n-\ — njj/n. Since the 
many body interaction is independent of the source of 
the polarization we expect our results to maintain their 
validity also in the case of an itinerant ferromagnet with 
a self-induced magnetic field, or when the polarization 
is achieved by other means, such as shining circularly 
polarized light on the sample. 

In section II, we study the large and the small wavevec- 
tor limits of the local field functions and their dependence 
with the electronic density and polarization. In section 
III, we give a simple parameterization of the local field 
factors which satisfy the asymptotic limits and repro- 
duces the most recent numerical results for the unpolar- 
ized electron gas.t3 Section IV presents our conclusions. 



II. LIMITING BEHAVIOR OF THE LOCAL 
FIELD FACTORS 

Our first approximation is to neglect the frequency de- 
pendence of the local field corrections. Althoj,Lfth the local 
field functions represent a dynamical effect, 
slowly on the scale of the Fermi frequency 
acceptable to neglect their frequency dependence if we 
are mainly interested in response functions. 
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A. Small wavevector 

At zero frequency and small wavevector, sum rules are 
used to connect the static limits of various response func- 
tions to certain thermodynamic coefficients, which can 
be expressed as derivatives of the ground state energy of 
the electron gasJl3 Subsequently, the renormalized local 
field functions, which are directly connected with these 
response functions, are written down as derivatives of the 
exchange and correlation energy of the interacting elec- 
tron gas (E xc ): 



2 sign{a' 



dr s d( 
(3) 
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2 dr s d( J ' 



where q = q/fcp is the normalized momentum, e xc = 
E xc /N is the exchange and correlation energy per parti- 
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FIG. 1: Initial slope a^(r a ,Q of the local field factor as a 
function of the polarization, (, for r 3 = 0.1 (dotted line), 
r s = 1. (dashed line), r a — 5. (dot-dashed line) and r s = 20. 
(solid line). 



cle measured in Rydbergs and the coupling strength, r s , 
is measured in units of the effective Bohr radius of the 
system.EJ Using the explicit expression of the exchange 
energy the local field functions become: 



G+(q^0) 



Gf(q-0) 



= ^(( 2 + 0v / T+C-Cv / W 
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(6) 



Similar expressions apply to Gf due to the fact that 

Gf(C) = G±(-C). 

The inclusion of the correlation energy in the small q 
limit of the local field factors is crucial to correctly eval- 
uate the response functions. If only the exchange contri- 
bution to the local field factors is included, the magnetic 
susceptibility of the-Lmpolarized gas develops a pole for 
/V2 ~ 2.221 Jill in disagreement with the results of 
extensive numerical calculations_that show a stable para- 
magnetic phase up to r s < 26.I1§Ej The addition of the 
correlation energy in the the small q limit of the local 
field factors prevents the occurrence of an instability in 
the unpolarized electron gas. 

To precisely evaluate the contribution from the corre- 
lation energy to G^(q — > 0) we use the latest numeri- 
cal computation of the ground-state energy of the twpj 
dimensional electron gas. Following Attaccalite et aL,Ej 



3 



0.7 



0.6 



-i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r~ 

s 

\ 

\ 

\ 

\ 



0.2 



0.5 - 



8 




_l I I I I I I I I I I I I l_ 



10 



15 



20 



0.8 



0.6 



0.4 



0.2 




FIG. 2: Initial slope a^(r s ,£) of the local field correction as 
a function of the coupling strength r s for different values of 
the spin polarization: f = — 1 (dashed line), £ = —0.5 (dotted 
line), £ = 0. (solid line), £ = 0.5 (dot-dashed line) and £ = 1 
(long-dashed line). 

the correlation energy per particle is given by: 

e c (r s ,() = (e- 0r °-l)e£\r s X) 

+a {r s ) + a 1 (r s )( 2 + a 2 (r s )( 4 , (7) 

where ei is the Taylor expansion of the exchange en- 
ergy beyond fourth order in £, (3 is a parameter and 
the functions. ai(r s ) are a generalization of previous 
expressions. E2I By taking the appropriate derivatives of 
Eq. (0), we derive analytical expressions for the ini- 
tial slope of the static local field correction functions: 

The behavior of c^f{r s ,Q and ct+(r s ,£) is quite 
different. Fig. [l] shows af(r s , £) as a function of the po- 
larization for several values of r s . Note that at is always 
positive, as the local effects always decrease the uniform 
electron density at large inter-electronic distances. Also 
note that for small values of r s is a monotonically in- 
creasing function of (. This behavior can be understood 
analyzing the contribution from the exchange interaction 
given by a+ = [(2 + C)\/TTC ~ C\/T Tr C]/27r. Since the 
exchange interaction takes place only between electrons 
with parallel spin, increasing values of £ induce larger ef- 
fects and further reduction of the uniform electron den- 
sity at small q. Consequently, the value of increases. 
The exchange contribution a+, which does not depend 
on r s , dominates at large electronic densities. However, 
correlation effects become more important with increas- 
ing r s and they partly cancel the strong £ dependence of 
a+, as it can be seen on Fig. |l|. 

Fig. H displays (r s , () as a function of r s for different 
values of the spin polarization. It can be noticed that 



FIG. 3: Initial slope a^(r 3 X) of the local field factor as a 
function of the polarization for different values of r 3 as indi- 
cated in the legend. 
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FIG. 4: Initial slope (r s , Q as a function of r s for different 
values of the spin polarization. 



is a monotonic increasing function of r s for values 
of C < 0.5. This feature is also displayed on Fig. 0. 
In particular, as previously noticed,t3 in the unpolarized 
electron gas a^(r s , £ = 0) increases with r s . 

Fig. H displays the initial slope aj(r s , £) as a function 

of C for different values of r s . As in the case of , the 
main contribution for small values of r s came from the ex- 
change energy: a J = — ( 2 J— g + 7= = )■ Thcrc- 

fore, for large electronic densities diverges for a fully 
polarized system (£ = ±1), behavior explained by the 



4 



fact that in a fully polarized gas the magnetic suscepti- 
bility becomes zero. As r s increases, correlation effects 
quench by a factor of e~@L? the diverging contribution 
from the exchange energy.tj Thus, correlation effects be- 
come dominant at small densities and the strong depen- 
dence of a~ with C is completely washed out. Instead, 
at large values of r s , becomes a linear function of £. 

Fig. |] displays a7 (r s , £) as a function of r s for different 
values of the spin polarization. The divergence at r s = 
and C = ±1 is clearly displayed. Note that is a 
monotonic decreasing function of r s for any value of C- 



B. Large wave vector 

In the limit of large wavevector it is easier to derive 
independently expressions for Gj and for G™. For large 
frequency or large wavevector, an iterative methiid gen- 
erates the exact asymptotic expressions for Gj .c3 For a 
two-dimensional electron gasjlil 



G+(q = /3+(r s , C) = 1 - - m (0) 



Gf(q-oo)=0f(r.,C) = 5ffu(O) 



(8) 
(9) 



where <7fj.(0) is the spin resolved pair distribution func- 
tion at the origin. It has been shown that ^(D-Lis 
largely unaffected by the degree of spin polarization. c3 cB 
In our calculation we use the simple expression <7f j_ (0) = 

7- ^-rr^z tt^z^t: — rr^ i where only the parameter 

[1 + 0.6032^ + 0.07263 r 2 s } 2 y v 

r s appears.t3 With this choice, /3f(r s ,Q = f)f(r a ,Q = 

(3+{r s ) and /3f(r s ,C) = ^(^,0 = P~ (r a ). 

In order to calculate the large momentum behavior of 
G™ we rieed to approximate the modified polarization 
functionlaO 
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f 
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■q,cr' 



+qer' 



(10) 



where £k<x = £k + sign(<r)7*i3 = k 2 /2m + sign(cr)7*_B is 
the quasiparticle energy in the static magnetic field B, 
Hk,<r is the exact occupation numbers in the interacting 
electron gas, and A is the area of the system. Since at 
large q the particle number renormalization is the dom- 
inant effect, Eq. ( |10| ) neglects the renormalization of the 
quasiparticle effective mass. 

By making an asymptotic expansion of the interact- 
ing polarization fupction,E£l, the local factor G™ can be 
written down as,E3E3 



G"(q->oc) 



At n 



V2(l + C) ; 



:q = l4r s ,Qq (11) 



where q — q/fcp is the normalized momentum, kp a = 
kpcr/kp the normalized Fermi momentum of the spin a 
electronic population and At a = i CT — t® is the difference 



between the kinetic energy of the electrons with spin a 

1 V- k 2 , . , 

m the interacting system, t a — — > , and m the 

N ^-^ Ira 

k 

I p 

non-interacting gas, t°„ = — nl- , over the total 

_/V ^— ' 2m 

k 

number of electrons. This equation is valid when At a 
is measured in Rydbergs and r s in units of the effective 
Bohr radius £3 

The difference in kinetic energies can be related 
with the exchange and correlation energies and their 
derivatives-, using the viriaB3 and the magnetic virial 
theorems :M 



Atn- = - 



r s d(e x + e c ) 
~2 ^ 



(12) 



where e x — E x /N and e c = E c /N are, respectively, 
the average exchange and correlation energy per parti- 
cle; t x a = E x /N and t c a = E^/N are the average exchange 
and correlation energy of the electrons with spin a. 
The average exchange energy for any spin population 

A /o 

is well known: t x = — (1 + sign(cr)C) 3/2 (Ry). The 

spin dependent correlation energies are difficult to eval- 
uate or to extract from numerical calculations. There- 
fore, we have to rely on same approximate scheme to 
extract the spin dependent correlation energy from the 
available computations of the full correlation energy. The 

E c 



N 



total correlation energy per particle is 

?l - ^ - + <f| - - . Perdew and WangJ 2 !! suggested the 

following parameterization for the correlation energy of 
a polarized electron gas: 



I (r. ) O=e c (r.,0) + fc(r„C)/(C) 



(13) 



where h(r s ,£) is an even function of the polariza- 

n + c\ A / 2 + fi - rW 2 - 2 

tion and /(C) = 1 — for a two- 

n .2(v2-l) 
dimensional systemJia This function can be decomposed 
as /(C) = [/ T (0 (l + + A(C)(1 - C)]/2, where f a = 
[a/1 + sign(a)£ — l]/[\/2 — 1]. An estimate of the spin 
dependent correlation energies can be obtained as: 



1 



where 



^(r s ,C)=e c (r s ,0) + 



MO , 
/(C) ' 



= (r s ,C)-e c (r s ,0)]. (14) 



Using the parameterization of the correlation energy 
proposed in Ref. |l^, Eq. (Q), we obtain reasonable val- 
ues for e c a . Also, its dependence with the polarization 
is the expected one. At a fixed value of r s , e'j (ep is a 
monotonic increasing (decreasing) function of C- 

Fig. | shows the dependence of 7| (r s , C) , Eq. (11), with 
the polarization for several values of r s . While the initial 
slopes (a^ ) an d the constant terms in the large q limit 
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FIG. 5: The coefficient 71 (r s , as a function of C, for different 
values of r a as indicated in the legend. 
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FIG. 6: The coefficient 7t(r s , £) as a function of the coupling 
strength r s for different values of the polarization: ( = —0.5 
(dashed line), £ = —0.25 (solid line), £ = 0. (long-dashed 
line), £ = 0.25 (dot-dashed line) and £ = 1 (dotted line). 



(Z^) of the local factors are always positive, the coeffi- 
cient 7 CT can have any sign. Since "f a oc Ato-, a change 
in sign of 7 CT implies a change of sign of At a = t a — t°. 
Even though, for a given spin Ai CT can be negative, the 
difference in the total kinetic energies of the interacting 
and the non-interacting system, AT = N[At^ + AiJ, is 
positive for any value of £ and r s Ej 

The negative value of Aij for negative polarizations 
is easily understood if the variation with the polariza- 
tion of the carrier effective mass is considered. It is well 
known that in a polarized electron gas the mass of the mi- 
nority (majority) spin populations-increases (decreases) 
with the value of the polarization.E3 Therefore, for large 



enough values of £ the mass renormalization dominates, 
and the kinetic energy of the minority (majority) carri- 
ers is reduced (increased) with respect to its free electron 
value. Also, as C increases, the exchange effects between 
the minority (majority) spins are greatly reduced (in- 
creased). Due to the combination of these two effects 
the electrons from the majority spin population will tend 
to accumulate around the few minority spins increasing 
the local charge around them. This generates a negative 
value of G™(q — > 00) for the minority spins. The opposite 
is true for the majority spins. 

At finite values of the polarization, the kinetic energy 
of each spin population is dominated by the exchange 
t s dc x 

contribution, At a = — — — e x ~ (/r s , a behavior that 

z OT s 

is independent of the parameterization used to obtain the 
spin-dependent correlation energy, Eq. ( |l4| ) in our case. 
At large values of r s , the correlation effects become more 
important as it can be seen on Fig. ^. Finally, we point 
out that the large effect of the polarization in the value 
of 7t(fs, C) is mainly due to the denominator (1 + £) 2 of 
Eq. Q. 

Fig. U displays 7|(r s , Q as a function of r s for several 
values of C between £ = —0.5 and £ = 1. The same trends 
as in Fig. || become apparent: strong dependence on £ 
and weak dependence with r s . For values of C, < 1, 7j 
increases with r s , reaches a maximum and decreases af- 
terwards. For example, for £ = the maximum is reached 
at r s = 2.88. For C = 1, 7| monotonically increases with 



III. SPIN DEPENDENT LOCAL FIELD 
CORRECTIONS 



Numerous parametrized expressions af the local field 
factors for j4jbe|p-unpolarizcd three-O'EjoH and two- 
dimensionalEZlL3E3 electron gas has been suggested 
since the pioneering work of Hubbard.EEl All these 
parametrized schemes rely on finding a function of the 
electronic density that smoothly interpolates between the 
short and long wavelength limits of the local field factors. 
However, parametrized expressions of for the polar- 
ized gas have not received so much attention. Numerical 
calculations in the polarized electron gas are also com- 
putationally challenging and, so far, only results for the 
ground state r £neijgy I -and the pair-distribution functions 
are available .EaOOH 

As we have discussed in the previous section, at small 
values of q the local factors follow a linear dependence, 
G^(q 0) = af(r s ,C) q, where q = q/fcp. In the 
opposite end of the spectrum the local corrections fol- 
low a linear plus constant dependence, G^(q — > 00) = 
/3 ± (r s ) -I- 7 CT (r s ,C) Q- Given the diversity of behaviors 
displayed by the parameters a^, ^ and 7^, we will con- 
sider a general interpolating scheme for all values of r s 
and We parametrized G ± with two fitting parameters, 
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TABLE I: Optimal fit 
parametrized in Eqs. 



parameters for the local field factors as 
@ and @. 
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FIG. 7: Local field correction factor G + (q, u> = 0) versus 
normalized momentum q/kF for r s = 1,2,5 and 10. Black 
circles correspond to the diffusion Monte Carlo (DMC) results 
of Ref. Ei| and solid lines are calculated according to Eqs. ( |li| ) 
and {ij^ 



q^ 1 and as: 



<£(q) 



afqe 



-Wit) 



(15) 



The fast decreasing exponential factors exp [— {q/q^') i ] 
and exp [— (g/c^) 4 ] are needed to mimic the rapid evo- 
lution of the local functions from their small q to their 
large q limiting behaviors. This fast evolution is clearly 
displayed on the latest diffusion Monte Carlo (DMC) 
results JL3 which show how the local factors G + (q) and 
G~(q) of the unpolarized gas follow the predicted lin- 
ear dependence at small values of the wavevector and 
rapidly reach their asymptotic large wavevector limits at 
not very large values of q. Since there are only numeri- 
cal results for the unpolarized electron gas, we consider 
the parameters q^ 1 and qf 1 as functions only on the elec- 
tron density. By fitting Eq. ( |i~5j ) to the DMC results 
of Ref. |l2| we find that the two fitting parameters are 
smoothly varying functions of the coupling strength and 
can be parametrized as: 



bfr. 



, ± 3/2 



± 3/2 



(16) 
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where the parameters af , bf and 4 are given in Table |. 



FIG. 8: Local field correction factor G~(q, u = 0) versus 
normalized momentum q/fc_F for the same values of r s . Black 
circles correspond to the diffusion Monte Carlo (DMC) results 
of RehjL^I and solid lines are calculated according to Eqs. ( |lE| ) 
and (16^ 



In Fig. [?] we compared our results for G + (q) of the 
unpolarized electron gas at r s =p-J , 2, 5 and 10, Eqs. jl5| ) 
and |[|), with the DMC results.Elln Fig.| our results for 
G~ (q) of the unpolarized electron gas are also compared 
with the numerical equivalent data. It is clear that our 
parametrized expressionS|-agree very well with the Quan- 
tum Monte Carlo results. E 2 ] By adding more free parame- 
ters to Eq. (|l|) the agreement with the numerical results 
will improve, but the question of how the parameters of 
the fit evolve with the spin polarization remains, as yet, 
unanswered. Thus, we keep the parameterization as sim- 
ple as possible: our fitting scheme, Eq. (jl5|), is the same 
for both local functions, it contains only two free param- 
eters and these two parameters have the same functional 
dependence with the electronic density, Eq. jlfj). 

Fig. ^ shows the momentum dependence of the local 
field functions Gj(q) for three values of the polariza- 
tion, ( = 0,0.5 and 0.9, and two values of the coupling 
strength, r s = 2 (density of 7.43 -10 10 cm~ 2 in GaAs) 
and r s = 10 (density of 0.297 TO 10 cm" 2 in GaAs). The 
factor associated with the majority spins, G^, behaves 
quite different from the one of the minority spins, Gf . 
The field factor G^ is always positive for positive values 
of the polarization. It slightly increases with the degree 
of polarization, but it keeps the characteristic peak of 
the unpolarized local factor around the same value of q. 
This peak-is a residue of the sharp peak in the exchange 
potential,EJ which is washed out and/or shifted to highp. 
values of q by the inclusion of short-range correlations.Ej 
On the other hand, the behavior of G^ is dominated by 
the forced change in its slope, from a positive value at 
small q to a large negative value at large wavevectors. As 
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+ 

(3 




degree of polarization. This is mainly due to the nega- 
tive value of the linear term on the large q limit of the 



FIG. 9: Local field corrections G^(q, u = 0) versus normal- 
ized momentum q/fcf for a two dimensional electron gas with 
r s — 2 (left panel) and r a — 10 (right panel). Results for 
G + at C, — (solid lines), Gj~ (long-dashed lines) and G|~ 
(dot-dashed lines) at £ = 0.5, and (dotted lines) and Gf 
(dashed lines) at £ = 0.9 are displayed. 



a result, G+ has always a maximum, whose position shifts 
to lower values of q with increasing £ and it is strongly 
dependent on the precise parameterization used. For ex- 
ample, it can shift to larger values of q if g^" and q^ 
are also functions of Q. This issue should be explored in 
greater detail. Finally, by comparing the field factors for 
r s = 2 and r s = 10 we conclude that their dependence 
on r s is weak. 

Fig. [Io]displays the local field factor G~ (q) versus nor- 
malized momentum for the same values of £ and r s used 
in Fig. ||. The main difference between G|~(q) and GJ (q) 
is that the latest one displayed a sharper peak around 
q ~ 2. It appears that higher order effects, which are 
important in the computation of Gt, cancel out in cal- 
culations of G7 due to the antisymmetric averaging over 
spin.E-3 It is also noticeable how the local factors change 
with increasing r s . 



IV. CONCLUSIONS 

We have considered an analytic parameterization of the 
spin dependent local field factors of the polarized two- 
dimensional electron gas, Eqs. ( |i~5| ) and Our pa- 
rameterization incorporates the known asymptotic limits 
of the local corrections and gives an accurate fit of the 
available quantum Monte Carlo data.E3 We found that 
the local field corrections associated to the minority spins 
strongly depend on the polarization, while the local field 
functions of the majority spins are less affected by the 
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FIG. 10: Local field corrections G^(q, u> = 0) versus nor- 
malized momentum q/fcf for a two dimensional electron gas 
with r s = 2 (left panel) and r 3 = 10 (right panel). Results 
for G~ at f = (solid lines), G^ (long-dashed lines) and GJ 
(dot-dashed lines) at £ = 0.5, and G^T (dotted lines) and GJ" 
(dashed lines) at £ = 0.9 are included. 



local field factors of the minority spins. 

The analytic parameterization used has only two pa- 
rameters which have been fitted to reproduce the latest 
DMC results for the unpolarized electron gas.liHl Since, 
that there are not data on polarized systems we have 
considered these two parameters as functions only on the 
electron density. Therefore, further study will be needed 
to evaluate their dependence on the spin polarization and 
the efficacy of our parameterization at large values of 

In conclusion, we believe that our approach provides 
a realistic qualitative description of the paramagnetic 
phase of the polarized electron gas. Caution, however, 
should be exercised in applying our calculation in the 
limit of C approaching unity, where the paramagnetic 
model breaks down. We have found that for small values 
of r s and large values of the magnetic susceptibility 
and the inverse dielectric constant develop a pole at the 
same value of the electronic density and the spin po- 
larization. This fact signals a charge-spin density wave 
instability in the polarized electron gas and it will be 
discussed elsewhereJlJ 
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